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Abstract 

The new simple method of the energy reconstruction for a com- 
bined calorimeter, which we called the e/h method, is suggested. 
It uses only the known e/h ratios and the electron calibration con- 
stants and does not require the determination of any parameters by 
a minimization technique. The method has been tested on the basis 
of the 1996 test beam data of the ATLAS barrel combined calorime- 
ter and demonstrated the correctness of the reconstruction of the 
mean values of energies. The obtained fractional energy resolution 
is [(58±3)%/\/^+(2.5±0.3)%]©(1.7±0.2)/^. This algorithm can 
be used for the fast energy reconstruction in the first level trigger. 



1 Introduction 



The key question of calorimetry generally and liadronic calorimetry 
in particular is the energy reconstruction. This question is especially im- 
portant when a hadronic calorimeter have a complex structure being a 
combined calorimeter. Such is the combined calorimeter of the ATLAS 
detector |^. 

In this paper we describe the new simple method of the energy recon- 
struction for a combined calorimeter, which we called the e/h method, 
and demonstrate its performance on the basis of the test beam data of the 
ATLAS combined prototype calorimeter. 



2 Combine Calorimeter 

The combined calorimeter prototype setup has been made consisting 
of the LAr electromagnetic calorimeter prototype inside the cryostat and 
downstream the Tile calorimeter prototype as shown in Fig. |I[ 
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The two calorimeters have been placed with their central axes at an 
angle to the beam of 12°. At this angle the two calorimeters have an active 
thickness of 10.3 A/. Beam quality and geometry were monitored with a 
set of beam wire chambers BCl, BC2, BC3 and trigger hodoscopes placed 
upstream of the LAr cryostat. To detect punchthrough particles and to 
measure the effect of longitudinal leakage a "muon wall" consisting of 10 
scintillator counters (each 2 cm thick) was located behind the calorimeters 
at a distance of about 1 metre. 

2.1 Electromagnetic Calorimeter 

The electromagnetic LAr calorimeter prototype consists of a stack of 
three azimuthal modules, each one spanning 9° in azimuth and extending 
over 2 m along the Z direction. The calorimeter structure is defined by 

2.2 mm thick steel-plated lead absorbers, folded to an accordion shape and 
separated by 3.8 mm gaps, filled with liquid argon. The signals are col- 
lected by Kapton electrodes located in the gaps. The calorimeter extends 
from an inner radius of 131.5 cm to an outer radius of 182.6 cm, repre- 
senting (at 7] = 0) a total of 25 radiation lengths (Xq), or 1.22 interaction 
lengths (A/) for protons. The calorimeter is longitudinally segmented into 
three compartments of 9 Xq, 9 Xq and 7 Xq, respectively. More details 
about this prototype can be found in [|l], 0. 

The cryostat has a cylindrical form with 2 m internal diameter, filled 
with liquid argon, and is made out of a 8 mm thick inner stainless-steel 
vessel, isolated by 30 cm of low-density foam (Rohacell), itself protected 
by a 1.2 mm thick aluminum outer wall. 

2.2 Hadronic Calorimeter 

The hadronic Tile calorimeter is a sampling device using steel as the 
absorber and scintillating tiles as the active material 0. The innovative 
feature of the design is the orientation of the tiles which are placed in planes 
perpendicular to the Z direction |Q. For a better sampling homogeneity 
the 3 mm thick scintillators are staggered in the radial direction. The tiles 
are separated along Z by 14 mm of steel, giving a steel/scintillator volume 
ratio of 4.7. Wavelength shifting fibres (WLS) running radially collect 
light from the tiles at both of their open edges. The hadron calorimeter 
prototype consists of an azimuthal stack of five modules. Each module 
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covers 27r/64 in azimuth and extends 1 m along the Z direction, such that 
the front face covers 100 x 20 cm^. The radial depth, from an inner radius 
of 200 cm to an outer radius of 380 cm, accounts for 8.9 A at ?7 = (80.5 
Xq). Read-out cells are defined by grouping together a bundle of fibres 
into one photomultiplier (PMT). Each of the 100 cells is read out by two 
PMTs and is fully projective in azimuth (with = 27r/64 ^ 0.1), while 
the segmentation along the Z axis is made by grouping fibres into read-out 
cells spanning AZ = 20 cm {Arj ^ 0.1) and is therefore not projective 
Each module is read out in four longitudinal segments (corresponding to 
about 1.5, 2, 2.5 and 3 A/ at = 0). More details of this prototype can be 
found in I, I, I, 

2.3 Data Selection 

Data were taken on the H8 beam of the CERN SPS, with pion and 
electron beams of 20, 40, 50, 80, 100, 150 and 300 GeV. 

We applied some similar to cuts to eliminate the non-single 

track pion events, the beam halo, the events with an interaction before 
LAr calorimeter, the electron and muon events. The set of cuts is the 
following: 

• the single-track pion events were selected by requiring the pulse 
height of the beam scintillation counters and the energy released in 
the presampler of the electromagnetic calorimeter to be compatible 
with that for a single particle; 

• the beam halo events were removed with appropriate cuts on the 
horizontal and vertical positions of the incoming track impact point 
and the space angle with respect to the beam axis as measured with 
the beam chambers; 

• the electron events were removed by the requirement that the energy 
deposited in the LAr calorimeter is less than 90 % of the beam energy; 

• a cut on the total energy rejects incoming muon; 

• to select the events with the hadronic shower origins in the first 
sampling of the LAr calorimeter; events with the energy depositions 
in this sampling compatible with that of a single minimum ionization 
particle were rejected; 
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• to select the events with the well developed hadronic showers energy 
depositions were required to be more than 10 % of the beam energy 
in the electromagnetic calorimeter and less than 70 % in the hadronic 
calorimeter. 



3 Existing Energy Reconstruction Methods 

Before going into description of the e/h method for the energy recon- 
struction we will first briefly review the existing algorithms used earlier for 
the energy reconstruction of the ATLAS combined prototype calorimeter. 

Three different algorithms have been developed in order to reconstruct 
the hadron energy of the ATLAS combined prototype calorimeter [0, llO]: 



the benchmark method 11 



the sampling weighting method, 



the cells weighting method ||T2| 



The benchmark algorithm is designed to be simple. With this method 
the incident energy is reconstructed with a minimal number of parameters 
(all energy independent with the exception of one). The energy of the 
particle is obtained as the sum of four terms: 

1. The sum of the signals in the electromagnetic calorimeter, EiAr, 
expressed in GeV using the calibration from electrons. 

2. A term proportional to the charge deposited in the hadronic calorime- 
ter, Rrue- 

3. A term to account for the energy lost in the cryostat, Edm ~ 

\ELAr,3 ■ O, ■ RTile,l\- 



4. A negative correction term, proportional to E\j^^. For showers that 
begin in the EM calorimeter, this term crudely accounts for its non- 
compensating behaviour. 

The reconstructed energy is: 

E = EiAr + a ■ RrUe + b ■ \/\ELAr,3 ■ 0- ■ RTile,l\ + C ■ El^r ■ (l) 
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The parameters a, b and c were determined by minimizing the fractional 
energy resolution a/E of 300 GeV pions. However the reconstructed energy 
is systematically underestimated. For this reason an additional step of 
rescaling is necessary. 

To rescale the reconstructed energy E to the beam energy Eheam the 
following expression was used: 

E, = , y ^ P , ^ , (2) 



F,^r[l~h{l-{^)^J]+Frue[l-h{l 



Tile 



where fh = I — f-j^o, f-^o is the energy-dependent fraction of the incident 
hadron energy, which is transferred to the electromagnetic sector, and given 



as a function of beam energy as in [^], two terms in the denominator are 



weighted by the average fractions of energy deposited in the LAr {Flat) 
and Tile {Fxne = 1 — -^LAr) calorimeters, taken from the data, the values 
of (|^)LAr and {^)Tiie (different for the two calorimeters) were found by 
fitting for all beam energies and turned out to be (^)lat- = 1-89 ± 0.04 



and (f^)riie = 1.22 ±0.02. As remarked in |]11|, the eh and eg values in the 



expression do not have the usual meaning as the calorimeter response 
to hadrons and electrons because they are determined not from the raw 
signals but from E, which already includes corrections. The rescaled mean 
values Er reasonably agreed with the nominal beam energies. 

The sampling weighting algorithm is based on a separate correction pa- 
rameter for each longitudinal compartment of the two calorimeters. These 
parameters are independently optimized for each incident energy and are 
indeed found to be energy-dependent. The correction strategy chosen here 
is to adjust downwards the response of the read-out cells with a large sig- 
nal to compensate for the response to large EM energy clusters, typically 
due to 7r° production. A separate weighting parameter was introduced for 
each longitudinal sampling. The energy measured in each readout cell E^ 
is corrected according to the formula: 

Er = E,-{l- W,^) , (3) 

where Ej is the energy sum over all cells of sampling j and Wj is the (pos- 
itive) weight to be optimized for each sampling j. In total eight energy- 
dependent parameters must be determined: one for each of the seven sam- 
plings, plus an additional conversion factor / to convert the hadronic signal 
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from charge to energy. It turned out that the samphng weighting technique 
improves the energy resolution but does not improve the hnearity obtained 
with the benchmark approach. 

The cells weighting method relies on correcting upwards the response 
of cells with relatively small signals, in order to equalize their response to 
that of cells with large (typically electromagnetic) deposited energies. The 
total energy is reconstructed correcting the energy in each cell of either 
calorimeter by a factor (typically > 1) which is a function of the energy in 
each cell and of the beam energy. The reconstructed energy is expressed 
as: 

= ^ Wem{Ecell, Ebeam) ' Ecell + ^ Whad{EceU, Eheam) ' Ecell + Edm, (4) 
cells cells 

Wem and Whad are the weights to be determined. The total number of 
parameters is reduced to 8 (including the cryostat constant). By minimiz- 
ing the energy resolution (with the constraint that the mean reconstructed 
energy reproduces the nominal beam energy), a value for each Wi is found. 

In all these methods the parameters are determined by minimization. 
So, it is necessary to have sufficiently large and unbiased sample of events 
and to perform the sophisticated off-line data treatment. Besides, for the 
two first methods it is necessary to rescale the reconstructed energy in 
order to obtain the nominal beam energy values. 



4 e/h Method of Energy Reconstruction 

4.1 Algorithm 

The response, Rh of a calorimeter to a hadronic shower is the sum of 
the contributions from the electromagnetic, i^e, and hadronic, E^, parts of 
the incident energy |n, |T3| 

Einc = Ee + Eh . (5) 

Rh = e-E, + h-Eh = e-Ei„,- (/,o + (h/e) ■ (1 - /,o)) , (6) 

where e (h) is the energy independent coefficient of transformation of the 
electromagnetic (pure hadronic, low-energy hadronic activity) energy to 
response, /^o = E^/Einc is the fraction of electromagnetic energy. From 
this 

e 1 

Einc = — ■ - ■ Rh , (7) 

TT e 
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where 

e e/h 



1 + {e/h - 1)^0 ■ 

For a combined calorimeter the beam energy, Eheam, deposits into the 
LAr compartment, EiAri the Tilecal compartment, E-rue, and into the 
passive material between the LAr and Tile calorimeters, Edm, 

Ebeam = EiAr + ErUe + Edm ■ (9) 

Using the expressions (|^ - (|) the following equation for the energy- 
reconstruction has been derived: 

E = CLAr ■ (e/ Vr)LAr ■ RlAt + CTile " (c/ T^)Tile " RtUc + Edm , (10) 

where 



(-) 



e/h) LAr 
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LAr 1 + {{e/h)LAr - l)Uo,LAr 



Uo,LAr = kiAr " lu E , (12) 

/e\ ^ {e/h)Tile . . 

W^^e 1 + ((e//i)„,, - l)/^o,™e ' ^ ^ 

/7r0,™e = knie " In (-Eriie) , (14) 



Edm — Cdm ■ y ELAr,3 ' Exile,! , (15) 
ELAr,3 = CLAr ■ (e/n) LAr ■ RLAr,3 (16) 

is the energy released in the third depth of the electromagnetic calorimeter, 

EtUs,! = CTile ■ {e/'7l)Tile " Rxile,! (17) 

is the energy released in the first depth of the hadronic calorimeter. In 
order to use the equation (|iy) it is necessary to know the values of the 
following constants which have been taken from [|l^ and equal to: CLAr = 

1/eLAr = 1.1, CTile = l/cTile = 0.145 ± 0.002, {e/h)LAr = 1.77 ± 0.02, 
{c/h)Tile = 1.3 ± 0.03, kLAr = ^Tile = 0.11, Cdm = 0.31. 

Note that the above formulae have been represented earlier in our work 
[ 161 devoted to the determination of the {e/h) ratio of the LAr electromag- 
netic compartment of the ATLAS barrel combined prototype calorimeter. 



8 



4.2 Iteration Procedure 



For the energy reconstruction by the formula (p!OD it is necessary to 
know the {e/7i)TUe ratio and the reconstructed energy itself. Therefore, 
the iteration procedure has been developed. 

Two iteration cycles were made. The first one is devoted to the deter- 
mination of the {e/Tr)Tiie ratio. The expression (^) can be written as 

/e\ ^ {e/h)Tiie ^^g^j 

KT^/TUe 1 + {{e/h)TUe - 1) " krUe ■ In {cTile " {e/7T)Tile " Rxile) 

As the first approximation, the value of {e/n)Tiie is calculated using the 
following equation 



(^Y = {e/h)TUe . . 

^71^ Tile 1 + {{e/h)Tile - 1) " krUe " In {cTile " 1-13 ■ Rnie) ' ^ ' 



where in the right side of this equation we used {e/7i)Tiie = 1-13 corre- 
sponding to fnO^Tiie = 0.5. The obtained value {e/7i)Tiie in i/-iteration 
(z/ = 0, 1, . . .) are used in (z/+ l)-iteration in the right side of the equation 
(]T8|) and the iteration continues. The iteration process is stopped when 
the convergence criterion 

I ie/7rr^,l - (e/TT)^,, | /(e/vr)^,, < e (20) 

is satisfied. 

The second iteration cycle is the determination of energy. As the first 
approximation, the value of E is calculated using the following equation 
with the {e/'K)Tiie ratio obtained in the first iteration cycle 

Eq = CLAr ■ (e/7r)°^^ • RlAr + CTile " {e/n)Tile " Rxile + Edm , (21) 

where 

/e\0 ^ (pLAr- 

Kt^JlAv 1 + {{f^)LAr - l)kLAr In {cLArl-27RLAr + CTile{^)TileRTile + ^dm) 

(22) 

In the right side of this equation we used (e/Tr)^^^ = 1.27 corresponding 
to f-n<-\LAr = 0.5 = 0.11 In (100 GeV). The convergence criterion is 

I E,+i -E,\/E,<e. (23) 
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Table 1: The average number of iterations < Na > for the various beam 
energies needed to receive the given value of accuracy e. 



E 


< Nu > 


< Nu > 


< Nu > 


<Nu > 


(GeV) 


e = 0.1% 


e = 0.3% 


e = 0.5% 


e = 1.% 


10 


1.79 


1.62 


1.49 


1.18 


20 


1.63 


1.39 


1.15 


1.10 


40 


1.41 


0.91 


0.78 


0.72 


50 


1.28 


0.80 


0.74 


0.69 


80 


0.77 


0.68 


0.58 


0.29 


100 


0.56 


0.30 


0.18 


0.09 


150 


0.82 


0.70 


0.62 


0.39 


300 


1.01 


0.74 


0.69 


0.63 



Table gives the average number of iterations < Nu > for the various 
beam energies needed to receive the given value of accuracy e. As can be 
seen, it is sufficiently only the first approximation for achievement , on 
average, of convergence with an accuracy of e = 1% for energies 80-150 
GeV and it is necessary to perform one iteration for the energies at 10 - 
50, 300 GeV. 

We specially investigated the accuracy of the first approximation of 
energy, Eq. Fig. ^ shows the comparison between the energy linearities, the 
mean values of E/Ebeam, obtained using the iteration procedure with e = 
0.1% (black circles) and the first approximation of energy (open circles). 
Fig. ^ shows the comparison between the energy resolutions obtained using 
these two approaches. As can be seen, the compared values are consistent 
within errors. 

The suggested algorithm of the energy reconstruction can be used for 
the fast energy reconstruction in the first level trigger. 

4.3 Energy Spectra 

Fig. ^ and Fig. ^ show the pion energy spectra reconstructed with the 
e/h method (e = 0.1%). The mean and a values of these distributions are 
extracted with Gaussian fits over ±2cr range. The obtained mean values 
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Table 2: Mean reconstructed energy, energy resolution and fractional 
energy resolution for the various beam energies obtained using the iteration 
procedure with e = 0.1% 





E (GeV) 


a (GeV) 


a/E (%) 


10* GeV 


9.30 ±0.07 


2.53 ±0.05 


27.20 ±0.58 


20* GeV 


19.44 ±0.06 


3.41 ±0.06 


17.54 ±0.31 


40 GeV 


39.62 ±0.11 


5.06 ±0.08 


12.77 ±0.21 


50 GeV 


49.85 ±0.13 


5.69 ±0.13 


11.41 ±0.26 


80 GeV 


79.45 ±0.16 


7.14 ±0.14 


8.99 ±0.18 


100 GeV 


99.10 ±0.17 


8.40 ±0.16 


8.48 ±0.16 


150 GeV 


150.52 ±0.19 


11.20 ±0.18 


7.44 ±0.12 


300 GeV 


298.23 ±0.37 


17.59 ±0.33 


5.90 ±0.11 



*The measured value of the beam energy is 9.81 GeV . 
*The measured value of the beam energy is 19.8 GeV . 



E, the energy resolutions a, and the fractional energy resolutions cr/E are 
listed in Table ^ for the various beam energies. 

4.4 Energy Linearity 

Fig. ^ demonstrates the correctness of the mean energy reconstruction. 
The mean value of E/Eheam equals to (99.5 ±0.3)% and the spread is ±1% 
except for the point at 10 GeV. But, as noted in Q, at this point the 
result is strongly dependent on the effective capability to remove events 
with interactions in the dead material upstream and to deconvolve the 
real pion contribution from the muon contamination. Fig. |^ and Fig. |^ 
compare the linearity, E/Ebeam, as a function of the beam energy for the 
e/h method, for the cells weighting method and for the benchmark method. 
For the cells weighting method the linearity is good to about ±1%. The 
benchmark method (without rescaling) badly reconstructs the beam 
energy: the reconstructed energy is systematically underestimated of about 
30%, the linearity depends from the beam energy and the ratio E/Eheam 
increases by 15% from 20 to 300 GeV. 



11 



4.5 Energy Resolutions 



Fig. p shows the fractional energy resolutions {cr/E) as a function 
of 1 / \/Eheam obtained by these three methods. As can be seen, the en- 
ergy resolutions for the e/h method (black circles) are comparable with 
the benchmark method (crosses) and only of 30% worse than for the cells 
weighting method (open circles). A fit to the data points gives the frac- 
tional energy resolution for the e/h method obtained using the iteration 
procedure with e = 0.1% 



a 
E 



(58 ± 3)% 



+ (2.5 ±0.3)% 



for the e/h method using the first approximation 
"(56 ±3)% 



a 
E 



(2.7 ±0.3)% 



for the benchmark method of 
"(60 ±3)% 



a 
E 



+ 1.8 ±0.2 



and for the cells weighting method of 
"(42 ±2)% 



a 
E 



1.8 ±0.1 



(1.7 ±0.2) 


GeV 


E 




ion (^ 




(1.8 ±0.2) 


GeV 


E 




(2.0 ±0.1) 


GeV 


E 




(1.8 ±0.1) 


GeV 


E 



(24) 



(25) 



(26) 



(27) 



where the symbol © indicates a sum in quadrature. As can be seen, the 
sampling term is consistent within errors for the e/h method and the bench- 
mark method and is smaller by 1.5 times for the cells weighting method. 
The constant term is the same for the benchmark method and the cells 
weighting method and is larger by (0.7 ± 0.3)% for the e/h method. The 
noise term of about 1.8 GeV is the same for these three methods that 
reflects its origin as the electronic noise. As to the two approaches for the 
e/h method the fitted parameters coincides within errors. 



5 Conclusions 

The new simple method of the energy reconstruction for a combined 
calorimeter, the e/h method, is suggested. It uses only the known e/h 
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ratios and the electron calibration constants and does not require the pre- 
vious determination of any parameters by a minimization technique. The 
method has been tested on the basis of the 1996 test beam data of the AT- 
LAS combined prototype calorimeter and demonstrated the correctness of 
the reconstruction of the mean values of energies. The obtained fractional 
energy resolution is [(58 ± 3)%/v^ + (2.5 ± 0.3)] (1.7 ± 0.2)/^;. This 
algorithm can be used for the fast energy reconstruction in the trigger of 
the first level. 
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Figure 2: Energy linearity as a function of the beam energy for the 
e/h method obtained using the iteration procedure with e = 0.1% (black 
circles) and the first approximation (open circles). 
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Figure 3: The fractional energy resolutions obtained with the e/h method 
(black circles), and the first approximation (open circles). 
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Figure 4: The energy distributions for Eheam = 10, 40 GeV (left column, 
up to down) and Ei^eam — 20, 50 GeV (right column, up to down). 
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Figure 5: The energy distributions for Eheam = 80, 150 GeV (left column, 
up to down) and -Bbeam — 100, 300 GeV (right column, up to down). 
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Figure 6: Energy linearity as a function of the beam energy for the e 
method (black circles) and the cells weighting method (open circles). 
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Figure 7: Energy linearity as a function of the beam energy for the e 
method (black circles) and the benchmark method (open circles). 
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Figure 8: The energy resolutions obtained with the e/h method (black 
circles), the benchmark method (crosses) and the cells weighting method 
(open circles). 
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